Astronomy & Astrophysics manuscript no. 7578 


© ESO 2008 


February 1, 2008 





Molecular gas in nearby low-luminosity QSO host galaxies 

T. Bertram 1 , A. Eckart 1 , S. Fischer 1 , J. Zuther 1 , C. Straubmeier 1 , L. Wisotzki 2 , and M. Krips 3 



o 
o 

(N 
C 



43 

6 
to 



> 
in 

(N 

o 
i> 
o 



X 



1 I. Physikalisches Institut, Universitat zu Koln, Ziilpicher Str. 77, 50937 Koln, Germany 
e-mail: bertram@phl.uni-koeln.de 

2 Astrophysikalisches Institut Potsdam, An der Sternwarte 16, 14482 Potsdam 

3 Harvard-Smithsonian Center for Astrophysics, SMA project, 645 North Aohoku Place, Hilo, HI 96720 



Received 2 April 2007 / Accepted 14 May 2007 



ABSTRACT 



Aims. This paper addresses the global molecular gas properties of a representative sample of galaxies hosting low-luminosity quasis- 
tellar objects. An abundant supply of gas is necessary to fuel both the active galactic nucleus and any circum-nuclear starburst activity 
of QSOs. The connection between ultraluminous infrared galaxies and the host properties of QSOs is still subject to a controversial 
debate. Nearby low-luminosity QSOs are ideally suited to study the properties of their host galaxies because of their higher frequency 
of occurrence compared to high-luminosity QSOs in the same comoving volume and because of their small cosmological distance. 
Methods. We selected a sample of nearby low-luminosity QSO host galaxies that is free of infrared excess biases. All objects are 
drawn from the Hamburg-ESO survey for bright UV-excess QSOs, have 6>-30° and redshifts that do not exceed z=0.06. The IRAM 
30m telescope was used to measure the 12 CO(1-0) and 12 CO(2-l) transition in parallel. 

Results. 27 out of 39 galaxies in the sample have been detected. The molecular gas masses of the detected sources range from 
0.4 • 10 9 M G to 9.7 • 1O 9 M . The upper limits of the non-detected sources correspond to molecular gas masses between 0.3 • 10 9 M Q 
and 1.2 • 10 9 M G . We can confirm that the majority of galaxies hosting low-luminosity QSOs are rich in molecular gas. The properties 
of galaxies hosting brighter type I AGN and circumnuclear starformation regions differ from the properties of galaxies with fainter 
central regions. The overall supply of molecular gas and the spread of the linewidth distribution is larger. When comparing the far- 
infrared with the CO luminosities, the distribution can be separated into two different power-laws: one describing the lower activity 
Seyfert I population and the second describing the luminous QSO population. The separation in the Lpm/L' co behaviour may be ex- 
plainable with differing degrees of compactness of the emission regions. We provide a simple model to describe the two power-laws. 
The sample studied in this paper is located in a transition region between the two populations. 

Key words, galaxies: active - galaxies: ISM - quasars: general 



1. Introduction 

The investigation of the molecular gas content and its distribu- 
tion in the host galaxies of quasistellar objects (QSOs) is a key 
issue in the understanding of evolutionary sequences and en- 
vironments of active galactic nuclei (AGN). All proposed se- 
quences involve star formation and an abundant supply of mate- 
rial to fuel the starburst and the central engine. 

ISanderseTai]dl988l) were the first to discuss a connection 
between ultraluminous infrared galaxies (ULIGs) and QSOs, 
suggesting that both represent different stages in an evolutionary 
sequence that starts with the collision of gas-rich spiral galax- 
ies. In the course of the merger, circumnuclear starburst activity 
is triggered and a dominant part of available molecular gas is 
concentrated within the center of the merging galaxy. Several 
interferometric studi es on individual UL IGs, like the archetypal 
ULIG Arp 220 (e.g. lScoville et al.ll 199 lb . reveal a massive con- 
centration of molecular gas in the center (SSI kpc diameter). The 
bulk of H2 in these gas rich objects seems to be associated with 
this central gas accumulation. This massive concentration not 
only favors starburst activity but may also be required to drive 
the nuclear activity that becomes apparent in the QSO phase of 
the sequence. The more distant and more gas-rich submillimeter 
galaxies (SMGs) resemble scaled-up versions of the local ULIG 
population, also showing compact CO emission regions confined 
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to the center. iTacconi et al.l {2006) state a median diameter of 
<4 kpc for 8 SMGs. 

Co ntrary to the evolutionary model, iDownes & Solomon] 
dl998h question the need for an AGN to power the FIR emis- 
sion in ULIGs. They argue that extreme starbursts in circum- 
nuclear molecular disks or rings are fully accountable for the 
high FIR luminosity. Following their line of arguments implies 
that, although several ULIGs with AGN are known, AGN are not 
mandatory for the ULIG appearance and ULIGs are not neces- 
sarily the predecessor s of QSOs. 

A recent study of Dunl op et alj ((2003) on the morphologi- 
cal properties of QSO host galaxies discusses a decrease of the 
number density of disk dominated host galaxies ending in a com- 
plete depletion at My < -23 (assuming Ho = 50 km s _1 Mpc" 1 ). 
From their analysis of Hubble space telescope and near infrared 
(NIR) data they conclude that QSO hosts have properties very 
similar to red, quiesc ent elliptical galaxies. These objects usu- 
ally are gas depleted. iBarthell d2006l) argues for a more differen- 
tiated picture at least of intermediate luminosity QSO hosts and 
points to several indications of recent or ongoing star formation 
activity in these objecs. 

Little is know n about mo lecular gas in nearb y QSO hosts: 
lEvans et"aT. I d200lLl2006l) and lScoville etail d2003)" discuss small 
samples of QSO host galaxie s selected from the Palo mar-Green 
(PG) Bright Quasar Survey dSchmidt & Greenlll983l) . Few stud- 
ies on individual sources add to the incomplete picture of molec- 
ular gas in local QSO hosts. A common denominator seems to be 
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Fig. 1. Lco-z distribution of already existing interferometric 
imaging studies of AGN hosts / interacting galaxies (thick 
outlined symbols) and the distribution of the nearby low- 
luminosity QSO sample. The underlying light gray symbols 
represent a compilation of the most CO lumin ous objects 
((U)L IGs/SMGs/HzRGs/QSOs) measured so far dGreve et alj 
120051) for the corresponding redshifts. The objects in the nearby 
low-luminosity QSO sample have redshifts that allow to re- 
solve structures on the sub-kpc scale in the 12 CO(1-0) transition 
with current mm-interferometers. [See the online edition of the 
Journal for a color version of this figure.] 



the presence of large amounts of molecular gas in the majority of 
these objects, which militates against a quiescent nature of host 
galaxies. The distribution of molecular gas, the extent and den- 
sity of regions emitting CO line emission and similarities with 
or differences to local Seyfert I galaxies have not yet been stud- 
ied on a solid statistical base. More detailed multi-wavelength 
investigations of a larger sample of nearby QSOs are certainly 
beneficial in the controversial debate on the nature and history of 
QSOs. Especially the separation of the starburst and the AGN 
component in extragalactic objects from the faint contribution 
of the underlying host galaxy require exceptionally high spatial 
resolution and sensitivity. The nearby QSOs with z<0.1, there- 
fore, represent an important link between the cosmologically lo- 
cal, less luminous AGN and the high redshift, high luminosity 
QSOs (at z>0.5). Within this volume the accessibility of impor- 
tant structural information with state-of-the-art interferometers 
is feasible. In the millimeter domain, interferometers then pro- 
vide sub-kpc resolutions that allow to probe the sizes of poten- 
tially compact circum-nuclear molecular gas reservoirs similar 
to the ULIG case (cf. Fig. Q]). If the molecular gas is less con- 
fined to the center compared to the ULIG case it will be possi- 
ble to study the distribution and potential signs of interaction on 
larger scales. 

We have selected a sample of nearby UV excess QSOs or 
luminous Seyfert I galaxies, that are not only suitable for inter- 
ferometric imaging in the millimeter wavelength domain but also 
allow for detailed imaging and spectroscopy in the NIR. For 41 
members of the sample single-dish CO data was obtained (in one 
case also interferometric data). The sample is introduced in Sect. 



[2] together with a definition of the term "low-luminosity QSO" 
that is used throughout this paper. The results of the observa- 
tions follow in Sect. [4] The data allows to draw conclusions on 
the total molecular gas content of low-luminosity QSOs and can 
be related to existing visible and far infrared quantities, as pre- 
sented in Sect. [5] Ho=75 km s Mpc -1 and qo=0.5 are assumed 
throughout the paper. 

2. The nearby low-luminosity QSO sample 

The only selection criterion for the sample of nearby low- 
luminosity QSOs was their small cosmological distance: only 
objects with a redshift z<0.060 were chosen. This redshift limit 
is based on a NIR spectroscopic constraint: it ensures the ob- 
servability of the diagnostic CO(2-0) rotation vibrational band 
head absorption line, which is important for the stellar popula- 
tion analysis. This line is then still accessible in the K-band. 

The members were selected entirely from an extended cat- 
alog of sources found in th e Hamburg/ESO survey (HES). The 
HES (Wis otzkietaT]|2000l) is a wide angle survey for optically 
bright QSOs, with a well-defined flux limit of Bj % 17.3, vary- 
ing from field to field, and a redshift coverage of 0<z<3.2. QSO 
candidates were identified in digitized objective prism plate data 
by applying a color- and spectral feature b ased selection scheme. 
With t he exception of objects listed in the lVeron-Cettv & VeronI 
(1996) catalog, all candidates were subject to spectroscopic 
follow-up observations to confirm the object's identity and to 
reject false classifications. 

The application of a starlike morphological criterion in many 
other QSO samples like the PG Bright Quasar Survey or the 10k 
catalog of the 2dF QSO Redshift Survey results in a significant 
degree of incompleteness at the low redshift end. One of the 
main advantages of the HES over other surveys is the consid- 
eration of extended objects. The HES sample, therefore, shows 
a high volume density of luminous type I AGN also at low cos- 
mological distances. This circumstance makes the HES sample a 
valuable source of objects for a study of nearby QSOs. A total of 
99 objects within the volume 0.01 <z<0.06 were identified by the 
HES. Thirty-nine of these form the subsample that is discussed 
in this paper. The objects in the subsample all have c5>-30°. To 
avoid potential selection biases, the composition of the subsam- 
ple matches the full sample of nearby low-luminosity QSOs in 
terms of redshift distribution and percentage of IRAS detected 
sources. 

It is important to note that no luminosity discrimination be- 
tween QSOs and Seyfert I galaxies was applied by the HES 
- all luminous type I AGN showing broad emission lines 
(FWHM>1000 km s~') in their follow-up spectra were included 
in their catalog. This has a direct implication on the abso- 
lute brightness distribution of the subsample here referred to 
as "nearby low-luminosity QSO sample". Our sample clearly 
probes th e low luminosity t ail of the local quasar luminosity 
function (Kohler et ldl 1 1997b . All objects in the sample have 
absolute Bj magnitudes exceeding (i.e., dimmer than) the tra- 
ditional boundarjQ M% — 22 between higher luminosity QSOs 
and lower luminosity Seyfert I galaxies. This boundary has no 
astrophysical motivation, a s it was technologically induced at 
the time of its introduction (Schmidt & Greenlfl9 83). However, 
to respect the commonly used definition of the term "QSO", we 
explicitely use the term "low-luminosity QSO" throughout the 
paper for objects identified in QSO surveys that may be fainter 
than the traditional boundary magnitude. Not only the HES but 
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also the PG Bright Quasar Survey provide low-luminosity QSOs 
in their samples. 

The restriction to type I AGN in the HES naturally resulted 
in the exclusion of type II AGN in the nearby low-luminosity 
QSO sample. For this reason the comparison with literature data 
in Sect. [5] was also restricted to type I AGN. 

Several members of the nearby low-luminosity QSO sam- 
ple were already subject to studies in the NIR and mm wave- 
length domain. First NIR imaging and spectroscopic results of 
nearby low-luminos ity QSO sample members are presented in 
iFischer et al.1 d2006lK A series of observations carried out with 
SEST and BIMA Bertram et al l 12006b preceded the observa- 
tions presented in this paper. In these first runs a slightly differ- 
ently composed sample was scanned for millimetric CO emis- 
sion. The primary goal was to identify the CO brightest ob- 
jects for high resolution interferometric follow-up observations. 
The resulting SEST spectra and CO properties of two detected 
nearby HES objects, HE 0108-4743 and HE 221 1-3903, are in- 
cluded in Sections [4] and [5] However, due to the higher detec- 
tion limit, they were not included in the FIR unbiased subsam- 
ple discussed in Sect. ED The third detection, HE 1029-1831, 
was included. Follow-up observations with the Plateau de Bure 
Interferometer allowed to resolve the central region in this ex- 
empla ry case. A detailed analysis can be found in iRrips et all 
J2007I) . 

3. Observations and data reduction 

in this paper we report on sensitive observations carried out with 
the IRAM 30m telescope on Pico Veleta (Spain) in September 
2005, June 2006 and January 2007. 

Twenty-six out of 38 observed sources were detected (cf. 
Table [TJ. In most cases the detection limits were much lower 
than in the preceding observing runs. However, part of the data 
obtained in June 2006 suffered from mediocre weather condi- 
tions and had to be excluded from further analysis. Moreover, in 
several cases the observing frequencies had to be based only on 
the redshifts provided by the HES. Although nominally uncer- 
tain by 300 km s _1 , in some cases the Hamburg/ESO based sys- 
temic velocities deviated from the measured systemic velocity 
centroids by up to 750 km s' 1 (cf. Sect. I5.21 , A certain fraction 
of weak sources may not have been detected because of this de 
facto uncertainty. 

The A and B receivers were used together with the 1 MHz 
and 4 MHz filter banks to acquire 12 CO(1-0) and 12 CO(2-l) 
emission line data in parallel. The configuration of the backends 
allowed for a velocity coverage Av between 1300 km s' 1 and 
1400 km s~' at the redshifted CO lines. The compactness of the 
sources allowed to take advantage of the wobbling secondary 
and the resulting excellent baselines. In the three observing runs 
useful data with a total of ~65h on-source integration time was 
obtained. 

The data reduction was carried out using the IRAM software 
CLASS/GILDAS. Each polarization direction was analyzed in- 
dividually and both directions were averaged whenever appro- 
priate. In all cases it was sufficient to subtract linear baselines. 
To increase the signal-to-noise ratio, most of the spectra (cf. 
Fig. |2]i were hanning smoothed. In order to obtain main beam 
temperatures, the main beam efficiency B e g = 0.75 and the for- 
ward efficiency F e fr = 0.95 were applied to T* A at ~ 1 10 GHz 
and fieff = 0.54 and F eff = 0.91 to T* A at -220 GHz. The 
line intensities Ico represent main beam temperatures integrated 
over the full velocity range of the emission line. The errors 
A/co were determined by the geometric average of the line error 



A/l = crvres yMl and the baseline error A/b = crv^N^I yNs, 
where <x is the channel-to-channel rms noise, v res the spectral 
resolution, Nl the number of channels over which the line is dis- 
tributed and A^b the number of channels used for the baseline 
fit. 



4. Results 

The CO line luminosity L co of the 12 CO(1-0) tran sition can be 
expressed as dSolom on. Dow nes. & Radford! 1992al) : 



L' co = 23.5 £Is*b D 2 L I co (l+zr 3 



[Kkm 



s- 1 pc 2 l 



i.e., as a function of the velocity-integrated line intensity Ico m 
units of K km s , the luminosity distance L\ measured in Mpc, 
and the solid angle of the source convolved with the telescope 
beam £2s*b m sc l uare arcseconds. For sources with solid angles 
much smaller than the telescope beam, Os+b can be approxi- 
mated with Qfieam- The expression L' CQ , commonly used in radio 
astronomy, and the general expression of the line luminosity Leo 
are related by L C o = (Snkv^Jc 3 ) L' CQ dSolomon et alJll992bb . 

IK km s~' corresponds to ~ 1.9 • 10 29 erg s~' at the 115 GHz 
transition. 

Established correlations (lYoung et al.l 1 1995b between the 
25.0 B mag arcsec' 2 isophotal diameter D25 of galaxies and the 
size of their CO emission region were used to confirm that, 
with the exception of HE 0212-0059 and HE 0232-0900, the 
sources should all be smaller than the ~22"telescope beam re- 
sulting at ~ 1 15 GHz. For the two exceptional cases it is likely 
that the overall CO line luminosity is underestimated in the 
prese nt study. In the case of HE 0232-0900, already published 
data dHorellou et al.l [l995l see below) allowed to determine a 
factor ~ 1.7 by which the current data is too low. Because of the 
very similar optical extent, a similar factor can be assumed for 
HE 0212-0059. 

The integrated line intensity, the systemic velocity (repre- 
sented by the flux- weighted centroid of the CO(1-0) line), the 
linewidth (FWZI), the resolution of the smoothed spectrum as 
well as the derived CO luminosity of each detected source is 
shown in Table [2] Also shown are estimates for the molecu- 
lar gas masses. They were obtained by applying a M(H2) / L' co 
conversion factor a to the CO luminosity. The determination 
of the H2 content using L' CQ as tracer has to be handled with 
care. Commonly used is a va lue of a=4.8 M Q (K km s~' 
pc 2 ) -1 dSolomon & Barrett! 199 lb . derived from galactic molecu- 
lar cloud observations. The common approach of using a 'stan- 
dard' conversion factor derived from galactic observations has 
the shortcoming of disregardin g the depend ency of a on the met- 
alicity of the region of interest (Israel 1997b. For a sample of ul- 
tra luminous infrared galaxies Downes & Solomon! £1998) state 
a conversion factor that is 5 times lower. Within this paper a — 4 
M Q (K km s _1 pc 2 )' 1 is adopted to allow for d irect compari- 
son wi th the results o flEvans etail d200lL[2006l) or lScoville et al.1 
(120031) . 

A subset of the sources had been subject to molecular gas 
studies before. The published data confirms the results listed 
in Table [2] The 12 CO(1-0) emission of the extended source 
HE 0232-0900 was m apped with the IRAM 30m telescope by 
iHorellou et al.l d 1995b . The value for Iqo of the central posi- 
tion agrees within 18%. After correction for the different beam 
sizes, the 12 CO(1- 0) detections of HE 0227 -0913 with the NRO 
45m telescope by IVila-Vilaro et al.l d 1998b . of HE 1237-0504 
with SEST bv ICurran et al.l d2000b . and of HE 0433-1028 with 
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Fig. 2. 12 CO spectra of detected host galaxies. In two cases, the weather conditions prevented reliable measurements of the 
12 CO(2-l) transition. 
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Fig. 2. — continued 



6 



T. Bertram et al.: Molecular gas in nearby low-luminosity QSO host galaxies 



Table 1. Journal of the observations carried out with the IRAM 30m telescope in in September 2005, June 2006, and January 2007. 
Three previously detected sources are also listed, z was taken from the HES. t mt denotes the total on-source integration time for each 
receiver. 



Obj. 


alt. name 


R.A. (J2000) 


Dec. (J2000) 


z 


Obs. 


Ant 






[h 


m 


s] 


[° 


' 


"] 




run 


[min] 


HE 0021-1810 


VCV2001 J002339.2- 175355 


00 


23 


39.4 


-17 


53 


53 


0.053 


•07 


118 


HE 0021-1819 


NPM1G -18.0010 


00 


23 


55.3 


-18 


02 


50 


0.052 


'07 


140 


HE 0040- 11 05 


VIII Zw 36 


00 


42 


36.8 


-10 


49 


21 


0.041 


'06 


112 


HE 0045-2145 


IRAS 00452-2145 


00 


47 


41.3 


-21 


29 


27 


0.021 


'05 


45 


HE 0108-4743 


IRAS 01089-4743 


01 


11 


09.7 


-47 


27 


36 


0.029 1 


SEST 


76 


HE 01 14-0015 


SDSS J01 1703.58+000027.4 


01 


17 


03.6 


+00 


00 


27 


0.046 


'06 


98 


HE 01 19-01 18 


II Zw 1 


01 


21 


59.8 


-01 


02 


25 


0.054 


•05 


60 


HE 0150-0344 


IRAS 01505-0343 


01 


53 


01.4 


-03 


29 


24 


0.046 


'06 


70 


HE 0203-0031 


Mrk 1018 


02 


6 


15.9 


-00 


17 


29 


0.043 


•06 


56 


HE 0212-0059 


Mrk 590 


02 


14 


33.6 


-00 


46 


00 


0.027 


'05 


60 


HE 0224-2834 


AM 0224-283 


02 


26 


25.7 


-28 


20 


59 


0.060 


'06 


55 


HE 0227-0913 


Mrk 1044 


02 


30 


05.4 


-08 


59 


53 


0.017 


•06 


46 


HE 0232-0900 


NGC 985 


02 


34 


37.7 


-08 


47 


16 


0.043 


•06 


32 


HE 0253-1641 


NPM1G-16.0109 


02 


56 


02.6 


-16 


29 


16 


0.032 


•05 


120 


HE 0345+0056 


IRAS 03450+0055 


03 


47 


40.2 


+01 


05 


14 


0.029 


•05 


117 


HE 0351+0240 


VCV2001 J035409.4+024931 


03 


54 


09.4 


+02 


49 


30 


0.034 


•06 


74 


HE 0412-0803 


MS 04124-0802 


04 


14 


52.6 


-07 


55 


41 


0.038 


•06 


115 


HE 0429-0247 


RXS J043 16-0241 


04 


31 


37.0 


-02 


41 


25 


0.041 


'06 


61 


HE 0433-1028 


Mrk 618 


04 


36 


22.2 


-10 


22 


33 


0.033 


•05 


25 


HE 0853-0126 




08 


56 


17.8 


-01 


38 


07 


0.060 


•06 


135 


HE 0853+0102 




08 


55 


54.3 


+00 


51 


10 


0.052 


•06 


120 


HE 0934+0119 


Mrk 707 


09 


37 


01.0 


+01 


05 


43 


0.051 


•06 


72 


HE 0949-0122 


Mrk 1239 


09 


52 


18.9 


-01 


36 


44 


0.019 


•05 


110 


HE 1011-0403 


PG 1011-040 


10 


14 


20.6 


-04 


18 


41 


0.057 


'06 


122 


HE 1017-0305 


Mrk 1253 


10 


19 


32.9 


-03 


20 


15 


0.048 


'05 


274 


HE 1029-1831 


NPM1G-18.0348 


10 


31 


57.3 


-18 


46 


34 


0.040 


PdBI 




HE 1107-0813 




11 


09 


48.5 


-08 


30 


15 


0.057 


'07 


80 


HE 1108-2813 


VCV2001 Jl 11048.0-283004 


11 


10 


48.0 


-28 


30 


03 


0.023 


'05, '07 


41 


HE 1126-0407 


Mrk 1298, PG1 126-041 


11 


29 


16.6 


-04 


24 


08 


0.060 


'05, '07 


41 


HE 1237-0504 


NGC 4593 


12 


39 


39.4 


-05 


20 


40 


0.009 


'05 


35 


HE 1248-1356 




12 


51 


32.4 


-14 


13 


17 


0.015 


'06 


102 


HE 1310-1051 


PG 1310-108 


13 


13 


05.7 


-11 


07 


42 


0.034 


'07 


64 


HE 1330-1013 




13 


32 


39.1 


-10 


28 


53 


0.022 


•07 


59 


HE 1338-1423 




13 


41 


12.9 


-14 


38 


40 


0.041 


'07 


85 


HE 1353-1917 


VCV2001 J135636.7-193144 


13 


56 


36.7 


-19 


31 


44 


0.034 


'07 


97 


HE 1417-0909 




14 


20 


06.3 


-09 


23 


13 


0.044 


'07 


106 


HE 2128-0221 


6dFJ2 130499-0208 14 


21 


30 


49.9 


-02 


08 


15 


0.052 


'05, '07 


295 


HE 2211-3903 




22 


14 


42.0 


-38 


48 


24 


0.039 


SEST 


114 


HE 2222-0026 


SDSS J222435.29-001 103.8 


22 


24 


35.3 


-00 


11 


04 


0.058 


'07 


225 


HE 2233-0124 


VCV2001 J223541.9+013933 


22 


35 


41.9 


-01 


39 


33 


0.056 


'07 


224 


HE 2302-0857 


Mrk 926 


23 


04 


43.4 


-08 


41 


09 


0.046 


'07 


150 



1 The redshift for HE 0108-4743 was take n from I Veron-Cettv & Veronl i 19961) and not spectroscopically verified by the HES. Since it differs 
significantly from other visual redshift data dLawrence et aljl999h . the value was excluded from the redshift discussion in Sect. 15,21 



SEST by IStrong et all (|2004|) agree within 10% with the val- 
ues i_gresej3^e^Jnthi^_gaper. HE 1237-0504 was also observed 
by lMaiolino et all d 19971) with the NRAO 12m telescope. Their 
beam-size corrected values for this object, for HE 0212-0059, 
and for HE 0949-0122, however, exceed the values presented 
here by a factor 2-3. For HE 0212-0059, the difference may be 
explained with the extension of the galaxy and the different beam 
sizes. The other results, however, are in contrast with the good 
agreements of the 4 other independent measurements mentioned 
above, which provide a high level of confidence in the data pre- 
sented here. 



5. Discussion 

5.1. Detection rate and molecular gas mass 

About 60% of the mostly southern sources in the nearby low- 
luminosity QSO sample have a declination > -30° and can be 
observed with the IRAM 30m telescope. Of these, about 65% 
have been subject to our study. They show a redshift distribution 
similar to the full nearby low-luminosity QSO sample (Fig. |3J. 
Within the restricted redshift range considered in this sample, the 
number of objects per redshift bin is purely volume dependent - 
the number increases with increasing z. When plotted against 
redshift (Fig. |4j>, the particular position of the sample in between 
previously studied local Seyfert I and more luminous QSO host 
galaxies becomes apparent. The omission of a far infrared (FIR) 
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Table 2. Summary of 12 CO properties. Unless otherwise noted, the data was acquired with the IRAM 30m telescope. 
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CO(l-O) 12 CO(2-l) 



Obj. 


v (LSR) 


D L 




AvpwZI 


Vies 


LS 2 


L^/10 8 3 


M(H 2 ) 




Avpwzi 


Vrcs 




[km/s] 


[Mpc] 


[K km/s] 


[km/s] 


[km/s] 




[K km/s pc 2 ] 


10"M G ] 


[K km/s] 


[km/s] 


[km/s 


HE 0021-1819 


15954 


215.5 


0.39±0.04 


110 


11.0 


A 


2.1 


0.8 


0.68 + 0.1 


121 


11.0 


HE 0040-1 105 


12578 


169.4 


0.96±0.08 


238 


10.8 


n 


3.2 


1.3 


1.19+0.15 


238 


22.0 


HE 0045-2145 


6403 


85.8 


8.20±0.18 


133 


9.0 


A 


7.2 


2.9 


22.39+0.59 


181 


5.3 


HE 0108-4743 4 


7175 


96.2 


2.24±0.07 


118 


7.4 


n 


10.4 


4.2 


3.67 + 0.1 


160 


7.4 


HE 01 14-0015 


13682 


184.5 


0.92±0.13 


326 


21.8 


n 


3.6 


1.5 


2.24+0.38 


305 


43.5 


HE 0119-0118 


16412 


221.7 


2.38±0.09 


165 


5.5 


A 


13.5 


5.4 


4.70+0.17 


143 


5.5 


HE 0150-0344 


14329 


193.3 


0.89±0.15 


305 


21.8 


n 


3.8 


1.5 


7 






HE 0212-0059 


7921 


106.3 


3.51 + 0.11 


240 


5.3 


n 


4.7 s 


1.9 s 


2.42+0.12 


251 


5.3 


HE 0224-2834 


18150 


245.6 


0.71+0.11 5 


132 


21.1 


A 


4.9 


2.0 


< 0.46 


132 


22.1 


HE 0227-0913 


4914 


65.8 


1.86±0.17 


254 


10.6 


n 


1.0 


0.4 


2.96+0.31 


254 


21.1 


HE 0232-0900 


12886 


173.6 


6.92±0.34 


597 


10.9 


n 


24.3 9 


9.7" 


8.13+0.66 


586 


22.0 


HE 0253-1641 


9580 


128.7 


2.89±0.15 


282 


13.4 


n 


5.6 


2.3 


3.89+0.37 


301 


21.0 


HE 0433-1028 


10651 


143.3 


9.42±0.27 


252 


3.4 


A 


22.6 


9.0 


16.55+0.48 


199 


5.4 


HE 0853-0126 


17899 


242.1 


1.28+0.15 


331 


22.1 


n 


8.6 


3.4 


0.82+0.36 


309 


44.1 


HE 0949-0122 


5905 


79.1 


1.23+0.14 


276 


21.2 


n 


0.9 


0.4 


3.55+0.24 


318 


21.2 


HE 1011-0403 


17572 


237.6 


1.37+0.12 


187 


11.0 


n 


8.9 


3.5 


7 






HE 1017-0305 


14737 


198.9 


1.25+0.10 


458 


21.8 


n 


5.7 


2.3 


2.44+0.23 


458 


21.8 


HE 1029-1831 6 


12112 


163.1 


4.24 ± 0.2 


254 


3.4 


A 


13.1 


5.3 


5.05 + 0.2 


254 


3.4 


HE 1107-0813 


17481 


236.4 


0.82+0.10 


308 


22.0 


n 


5.2 


2.1 


<0.3 


310 


44.0 


HE 1108-2813 


7198 


96.5 


8.32+0.28 


309 


5.3 


n 


9.2 


3.7 


18.84+0.47 


320 


5.3 


HE 1126-0407 


18006 


243.6 


1.52+0.15 


640 


22.1 


n 


10.3 


4.1 


3.02+0.20 


551 


22.1 


HE 1237-0504 


2531 


33.8 


7.52+0.57 


302 


13.1 


n 


1.0 


0.4 


6.04+1.01 


315 


21.0 


HE 1248-1356 


4338 


58.0 


2.96+0.19 


348 


10.5 


n 


1.2 


0.5 


4.23+0.39 


338 


21.1 


HE 1330-1013 


6744 


90.4 


1.82+0.16 


266 


10.6 


n 


1.8 


0.7 


2.4 + 0.2 


245 


10.6 


HE 1353-1917 


10472 


140.8 


3.25+0.22 


592 


10.7 


n 


7.6 


3.0 


5.5 + 0.3 


549 


10.8 


HE2211-3903 4 


11945 


160.8 


1.87+0.11 


466 


30.0 


n 


23.2 


9.3 


1.85+0.06 


260 


15.1 


HE 2222-0026 


17414 


235.5 


0.21+0.04 


154 


22.0 


A 


1.3 


0.5 


0.61+0.06 


176 


22.0 


HE 2233+0124 


16913 


228.6 


0.89+0.08 


506 


22.0 


n 


5.4 


2.1 


0.95+0.12 


506 


22.0 


HE 2302-0857 


14120 


190.4 


2.31+0.16 


653 


10.9 


n 


9.7 


3.9 


5.17+0.29 


817 


10.9 


Non-detections 


HE 0021-1810 


16039 


216.6 


< 0.07 


315 


43.8 




< 1.1 


< 0.5 


< 0.12 




43.8 


HE 0203-0031 


12723 


171.4 


< 0.16 


315 


43.4 




< 1.6 


<0.6 


< 0.51 




43.4 


HE 0345+0056 


8994 


120.8 


< 0.15 


315 


42.9 




< 0.8 


<0.3 


<0.26 




42.9 


HE 0351+0240 


10793 


137.0 


< 0.19 


315 


43.1 




< 1.2 


<0.5 


7 






HE 0412-0803 


11392 


152.9 


< 0.14 


315 


43.2 




< 1.2 


< 0.5 


7 






HE 0429-0247 


12441 


165.5 


< 0.31 


315 


43.4 




< 2.9 


< 1.2 


7 






HE 0853+0102 


15589 


210.5 


< 0.16 


315 


43.8 




< 2.4 


< 1.0 


7 






HE 0934+0119 


15091 


203.7 


< 0.12 


315 


43.7 




< 1.7 


<0.7 


< 0.24 




43.7 


HE 1310-1051 


10193 


137.0 


< 0.16 


315 


43.0 




< 1.0 


< 0.4 


< 0.23 




43.0 


HE 1338-1423 


12528 


168.8 


< 0.15 


315 


43.4 




< 1.5 


<0.6 


< 0.21 




43.4 


HE 1417-0909 


13191 


177.8 


<0.09 


315 


43.4 




< 1.0 


< 0.4 


< 0.18 




43.4 


HE 2128-0221 


15828 


213.8 


< 0.08 


315 


43.8 




< 1.3 


< 0.5 


< 0.12 




43.8 



1 Upper limits and errors represent lcr. 2 Line shape. A denotes a triangular shaped line profile, as it is expected for turbulent line-of-sight velocity 
dispersion, n represents boxy or double-horn line profiles indicative for emission from an inclined, rotating disk. 3 For the upper L' co limits, the 
average linewidth Av = 315km s and 3cr detection limits were considered. 4 SEST 15m telescope measurement. A beam efficiency B e ff = 0.7 
was applied to and a 45" beam was considered for L™. 5 3cr detection level for T mb . 6 IRAM PdBI measurement (Krips et al. 2007), resulting 
inSAv = (21 ± l)Jykms-' for 12 CO(1-0) and 5 Av = (25 ± l)Jy km s~' for 12 CO(2-l). S/T mh = 4.95 Jy/K was used to compare the PdBI results 
with the 30m telescope data. 7 1mm data not reliable due to high water vapor. 8 Source is likely to be extended. Value probably too low by an 
estimated factor 1.7 9 Source is not compact. Values too low by a factor ~1.7 



selection and the on average larger distance compared to the lo- 
cal Seyfert I population resulted in the inclusion of several ob- 
jects with FIR flux densities below the IRAS detection limit (cf. 
Sect. 15.41 ). In several of these cases we still were able to detect 
CO emission. 

Our 3cr detecti on limits agree with the lim its in the PG 
QSO h ost studies bv lEvans et all d200lLl2006h and lScoville et alj 
(2003), carried out with the OVRO array and the IRAM 30m 



telescope. Differ ent from Evans et al. d200lll2006h . the volume 
limited sample of lScoville et alj d2003l) . consisting of 12 objects, 
is not confined to FIR selected sources. Therefore, a direct com- 
parison with their results is appropriate. At redshifts exceeding 
0.06, the common detection limit prevents the detection of galax- 
ies with a molecular gas content of below 10 9 M - only gas 
rich objects are then detectable. The redshifts of v irtually all PG 
QSOs exceed this value. Nevertheless, sample of lScoville et alj 
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Fig. 3. The number of existing / observable with the IRAM 30m 
telescope / observed / detected HE QSOs with increasing vol- 
ume. At the distances discussed here, the frequency of occur- 
rence is dominated by the size of the comoving volume. 



10" 



10" 



10° 



nearby HE QSOs (this sample) 
PG QSOs (Evans et al. 2006) 
PG QSOs (Scoville et al. 2002) 
Seyfert I Galaxies (Heckman etal.~T889) 
Seyfert I Galaxies (Maiolino et al J997| 
Seyfert I Galaxies (SJtrong et al. 2004) " " 
LIGs (Sanders et al. 19 + 91) + +* 
ULIGs (Solomon etal. 1997) ++ + + 




0.001 



Fig. 4. CO luminosity versus redshift. For objects represented 
by filled shapes Lfir is available, grey shapes indicate available 
upper limits for L FIR and outlined shapes represent objects that 
have not been listed in the IRAS FSC. The shaded area in the 
bottom right corner indicates the region in which the determi- 
nation of Ico with the 30m telescope falls below 0.315 K km 
s _1 . This value corresponds to a signal of 1 mK with a linewidth 
of 315 km s , the average width of the detected lines. [See the 
online edition of the Journal for a color version of this figure.] 



d2003l) resulted in a detection rate of 75% which lets the authors 
conclude that the majority of luminous, low-redshift QSOs have 
gas-rich host galaxies. This picture may become more compli- 
cated when considering the incompleteness of the PG survey at 
low redshifts. The rejection of extended objects may cause an 
underrepresentation of bulge-dominated galaxies with a lower 
gas content. On the other hand, due to the arbitrariness of the 

Mg 22 borderline, low-luminosity QSOs above the line still 

may have host properties that are similar to their brighter equiva- 
lent. With 70% the detection rate of the IRAM 30 m observations 
presen ted here is almost identical to the rate of IScoville et alj 
(2003). The molecular gas masses of the detected host galaxies 




-900 -750 -600 -450 -300 -150 150 300 450 
Velocity offset [km s" 1 ] 

Fig. 5. The distribution of the differences between the CO veloc- 
ity and the velocity derived from restframe visible spectroscopy. 
The systematic blueshift of the spectroscopic features deter- 
mined in the visible is an indication for outflows or partial dust 
obscuration in the vicinity of the active nucleus. 



range from 0.4T0 9 M o to 9.7-lO 9 M with an average molecular 
gas mass of 2.8TO 9 M (3.0T0 9 M o when including the 2 sources 
that were detected with SEST). 

5.2. Comparison of CO and optical redshifts 

Millimetric CO observations are well suited for the accurate 
measurement of systemic velocities of active galaxies. CO traces 
the cold and extended molecular gas distribution of the total host 
galaxy rather than the reservoir of highly excited gas in proxim- 
ity to the central source. Features of the latter are commonly 
used to determine the redshift of active galaxies in the visible 
wavelength domain. These redshifts, however, can be affected 
by outflows, dust obscuration or other asymmetric phenomena 
that result in apparent velocities deviating from the actual sys- 
temic velocity. 

For nearby objects like the sample discussed in this paper, 
the redshifts determined in the course of the HES followup spec- 
troscopy progr am are mainly base d on the narrow J O ni|/t5007 
emission line dReimers et alJll99q) . iBorosonl d2005l) shows that 
the [Om] line is blue-shifted relative to the systemic velocity 
in up to 50% of all AGN. The emission peak shift may be as 
large as several hundred km s _I . Fig. |5]shows the distribution of 
offsets between the velocities corresponding to the HE redshifts 
and the CO velocities measured in this study. Negative velocity 
differences imply blue-shifted visible emission features. On top 
of the nominal uncertainty of HE velocities (~300 km s _1 ), the 
nearby low-luminosity QSOs indeed show a tendency for such a 
blue-shift with a mean velocity offset of -174 km s _1 . 

5.3. Blue luminosity and CO line characteristics 

Fig. [6] shows the connection between absolute Bj brightness 
and CO luminosity. The Bj fluxes have a photometric uncer- 
tainty of ~0.2 mag and were determined by the HES group, 
based on their spectral plate data. Each value represents the flux 
within the central seeing disk, which for type I AGN at these red- 
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Fig. 6. Absolute B j magnitude versus CO luminosity. The mag- 
nitudes shown here represent the central seeing disk, includ- 
ing contributions of both AGN and star formation in the cen- 
tral part of each galaxy. The majority of objects showing M(B j) 
brighter than -20 mag are also CO luminous. Exceptions are (1) 
HE 0345+0056, (2) HE 1338-1423, and (3) HE 0203-0031. 
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1 Objects with M(Bj) > -20 mag 
3 Objects with M(Bj) < -20 mag 
- total 
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Fig. 7. Distribution of CO line widths (FWZI). Linewidths ex- 
ceeding 400 km s _1 can be commonly found only in bright ob- 
jects with absolute blue magnitudes < -20. 

Table 3. Mean CO(1-0) line widths and standard deviations of 
the mean for different line shapes and different absolute blue 
magnitudes. The numbers in parentheses represent the number 
of objects that fall into the corresponding category. 









A 




n 


a + n 


M(Bj) > 
M(Bj) < 
all 


-20 
-20 


132 
200 
171 


± 13 (3) 
± 28 (4) 
± 22 (7) 


301 
426 
364 


±34(11) 
±49(11) 
± 32 (22) 


265 ± 33 (14) 
367 ± 45 (15) 
317 ± 29 (29) 



shifts is dominated by the nuclear contribution, with only small 

proportions from the host galaxy. 

A more detailed description is given in Wis otzki et alj 
(2000). The magnitudes shown here are not corrected for ex- 
tinction and, therefore, systematically too faint by a mean value 
of ~0.23 mag. 

The data does not allow to distinguish between host and 
AGN contribution. Thus an increase of By luminosity can be 
the result of enhanced star formation activity in the central 1- 
2 kpc, enhanced AGN activity or both. In any way, the increase 
in activity comes along with larger reservoirs of molecular gas, 
as can be seen in Fig. [6] Except for 3 non-detected sources, low- 
luminosity QSOs with absolute central magnitudes brighter than 
M{Bj) = -20 mag all have CO luminosities of > 5 • 10 8 K km 
s _1 pc 2 , while the majority of fainter objects contain less molec- 
ular gas. This correlation seems to be intrinsic and not the result 
of a distance selection bias (with more luminous AGN being at 
larger distances with higher CO detection limits), as the small 
fraction of non-detections in the bright case shows. It remains 
unclear, whether the 3 By luminous non-detections are actually 
gas depleted or whether they were not detected because of an 
uncertain CO redshift. Upcoming HI observations, but also anal- 
yses of the morphologies in the visible/NIR will provide further 
information on the nature of these objects. 

Not only the total CO content is higher but also the spread 
of the CO(1-0) line widths is wider for the brighter AGN hosts. 
Fig.|7]shows the distribution of line widths in all detected galax- 
ies. With one exception, all objects showing FWZI line widths 
broader than 400 km s _1 have absolute central magnitudes of 
M(Bj) < -20 mag. The average line width of these bright ob- 
jects is 367 km s _1 , ~ 100 km s broader than the fainter sources 
(cf. Tab. |3]l. The average signal-to-noise ratio (S/N) of the 
CO(1-0) detections in the B j luminous case is only slightly (1.2 
times) larger than the average S/N in the fainter case. Therefore, 
a significant S/N bias that potentially favors the identification of 
broader line profiles in more active objects is not present. 

Such a different behavior could result from a selection effect 
in the optical identification process of QSO candidates. Objects 
with faint AGN activity may only be identified as type I AGN 
host when seen almost face-on, leaving the Seyfert I nucleus un- 
obscured, while brighter AGN contributions can also be identi- 
fied in objects with higher inclination. A scenario of an intrin- 
sically broader line width distribution, however, cannot be ruled 
out. A separation of sources into a class of triangular shaped and 
a second class of box shaped or double horn line profiles resulted 
in similar counts for each class (cf. Tables [2] and [3} in both the 
faint and the bright cases. The triangular shaped lines are indica- 
tive for a face-on view onto a turbulent gas distribution without 
a dominating disk component in the line-of-sight direction. The 
second class represents objects with rotating CO emission disks 
or rings at a certain inclination with respect to the line-of-sight. 
In several low-signal cases the assignment of a class was not ob- 
vious. Nevertheless, an aforementioned selection effect that fa- 
vors low-inclination views on faint AGN candidates should yield 
a higher fraction of triangular shaped line profiles, which is not 
the case. Noting the caveat of small number statistics, we would 
like to point out the difference in the mean line widths of trian- 
gular shaped line profiles. Bright objects with triangular shaped 
line profiles have a mean line width of 200 km s _1 , whereas in 
the faint cases the mean line width is only ~130 km s . This 
difference could be a sign of higher velocity dispersion in the 
molecular gas distribution of the more active objects. 

Most of the line profiles do not show strong asymmetries that 
can be used as indication for ongoing interaction. Clear asymme- 
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Table 4. Far infrared properties of the IRAS detected sources in 
the subsample. The 60 and 100 ym IRA S flux densities are taken 
from the IRAS Fai nt Source Catalog dMoshir & et al.ll 19901) . or 
in few cases from iRice et alj 0988), San ders et al.1 d 1989b . or 
ISanders et all d2003l) . 



Obj. 










[Jy] 


[Jy] 


[10 10 L o ] 


HE 0045-2145 


3.6 


5.3 


4.2 


HE 0108-4743 


1.0 


2.2 


1.7 


HE 0119-0118 


1.5 


1.8 


10.8 


HE 0150-0344 


0.5 


<1.7 


<4.4 


HE 0212-0059 


0.5 


1.5 


1.2 


HE 0224-2834 


0.4 


0.4 


3.3 


HE 0227-0913 


0.4 


0.9 


0.3 


HE 0232-0900 


1.4 


1.9 


6.5 


HE 0253-1641 


0.7 


0.8 


1.6 


HE 0345+0056 


0.5 


<3.2 


<2.7 


HE 0412-0803 


0.6 


<1.4 


<2.8 


HE 0433-1028 


2.7 


4.2 


9.1 


HE 0949-0122 


1.3 


<2.4 


<1.5 


HE 1011-0403 


0.2 


<0.3 
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tries can be seen in the two sources that are probably larger than 
the telescope beam (HE 0212-0059 and HE 0232-0900). In these 
cases pointing errors are the likely explanation for the asymmet- 
ric line profile. 
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Fig. 8. FIR luminosity as a function of CO luminosity for the 
nearby QSO sample and various other Seyfert I, PG QSO host, 
luminous infrared galaxies and normal galaxies. The solid line 
represents a power law fit to the PG QSO hosts, t he dashed lin e 
the distance limited sample of normal galaxies of ISagd (Q993), 
the dotted line a flux-limit ed sample of normal spiral galax- 
ies (ISolomon & Sagdll988l) . and the dash ed-dotted line repre - 
sents a power-law fit to the LIG sample of Sand ers et alj dl99ll) . 
Outlined symbols represent upper L' CQ limits. [See the online 
edition of the Journal for a color version of this figure.] 



5.4. Far infrared luminosity and star formation efficiency 

The correlation between molecular gas and far infrared emission 
is well known and often discussed in the context of CO stud- 
ies. Lfir/Xp or Lfir/M(H2) is commonly referred to as Star 
Formation Efficiency (SFE) indicator, assuming that L' CQ traces 
the cold, gravitationally bound molecular gas reservoirs that 
form the birth places of young, dust enshrouded stars. However, 
CO is only an indirect measure as it actually traces the overall 
molecular gas content in galaxies rather than explicitly the dens e 
cores that produce young stars (e.g., iGao & Solomonl [2004). 
Furthermore, the far IR emission of galaxies is a composite of 
an active star formation and a quiescent cirrus-like component 
dHeloulll986l) . In the case of AGN hosts, AGN driven dust heat- 
ing may also contribute significantly to the FIR luminosity. 

Many samples are based on IR selection criteria and, there- 
fore, represent only the IR and CO bright tail of the lu- 
minosity distribution. So do several stu dies of the molecu- 
lar gas content of Seyfert galaxies (e.g. Heckman et"ail l 1989; 
iPapadopoulos & Seaquisdll998tlStrong et al.ll2004l) or the study 
of PG QSO host galaxies by lEvans etal J d200lll20"06l) . For the 
sample discussed in this paper an IR selection criterion was not 
introduced. Only ~50% of the nearby HE QSOs sample are 
listed in IRAS catalogs. The subsample discussed here contains 
the same fraction of IRAS sources (cf. Table |4]l - an IR bias 
was carefully avoided. 65% of the detected sources are listed 
in IRAS catalogs, while only 2 of the 12 non-detections are de- 
tected by IRAS. Only HE 2302-0857 is located in a field that 
was not scanned by IRAS. 



Many of the IRAS detected sources only show 60 //m and 
100 yum flux densities and several of them are close to the IRAS 
detection limit. Due to the limited data at hand, the following 
discussion will focus on Lfir. It is based on the 60 fim and 
100 jum IRAS flux densities and represents the luminosity 
between 43/vm and 123/rni. Due to this bandwidth limit, warmer 
components peaking at shorter wavelengths, especially potential 
AGN heated dust components remain less considered in Lfir. 
In the context of LIGs/ULIGs or nearby galaxies, Lir is often 
used instead of Lfir- Lir also includes the 12//m and 25^m 
IRAS bands and represents the total IR luminosity ranging 
from 8 to 1000 yum. Fo r the definition of Lfir and Lir cf. 
ISanders & Mirabelld 19961) . 

Fig. |8] shows the distribution of the IRAS detected sources 
in the L' CQ - Lfir diagram, together with the C O data of a 
volume limited sa mple of nearby spiral galaxies (|Sagd[l 993). 
Seyfert I galaxies dHeckman et al.1 119891: iMaiolino et al.i l 1997; 
IStrong et al.ll2004l) , luminous (LIGs. ISanders et al.l 199 iFan d ul- 
traluminou s infrared galaxies (U LIGs, Solo mon et al.l 1997b . and 
PG QSOs (IScoville et alj|2003t lEvans et al.ll2006l) . Also 'shown 
are linear regression fits (in log-log space) to normal and to lu- 
minous IR galaxy samples as well as to the PG QSO hosts. 

The distribution of type I AGN hosts in the plot indicates the 
presence of two populations with differing power-laws: one pop- 
ulation that follows the power-law of normal spiral galaxies and 
a second population that has Lfir/L^ properties very similar to 
LIGs. 
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The studies of lSolomon & Sagd (11988b and|Sage| d 1993b are 
used to quantify the power law for normal galaxies: A fit to the 
distance limited sa mple of norm al galaxies (excluding the upper 
limits considered in Sagel dl993l) ) resulted in 

L FIR /10 9 = 8-[L co /10 9 ] a72(±005) , 

while a fit to a FIR-selected sample of isolated spiral galax- 
ies or spiral galaxies t hat do not show signs of interaction 
dSolomon& Sagdl 19881) yielded 

L FIR /10 9 = 31 • [L co /10 9 ] a74(±007) . 

Due to the FIR selection criterion, the latter sample is biased 
towards more active star formation with respect to average spiral 
galaxies. The bulk of local Seyfert I galaxies follows the trend 
described by these two fits, wi th a somewh at higher Lfir/L co 
ratio for any given L' CQ than the Sage ( 1993) sample. 

The power law fit to the PG QSO hosts, on the other hand, 
produces a slope of about unity: 

L FIR /10 9 = 75-[L CO /10T 96(±OI9) 

and is almost identical to the fit to the LIG data: 

iFiR/10 9 = 71-[L^ o /10 9 ] - 96(±009) . 

This common trend may be indicative for similar star formation 
activity in LIGs and luminous QSO hosts. 

As for the low-luminosity QSO sample, Figure [8] supports 
the picture of the sample being a link between the local Seyfert I 
population and luminous QSO host galaxies. The sample mem- 
bers seem to follow both trends, most of them can be found in a 
transition region between the two populations. 

5.5. Opaque CO and FIR emission in (U)LIGs and QSO host 
galaxies? 

In the following we will summarize the various published re- 
sults on the L' CQ -Lpi R relation in normal and infrared luminous 
galaxies together with the findings on the nature of the emission 
regions. We combine different pieces of information to a simple 
qualitative model that explains geometrically the separation of 
type I AGN hosts into two populations in the Lfir range between 
1O 9 L and 5 ■ 10 n L G . 

As discussed in Sect. 15 .41 normal, undisturbed galaxies show 
a nonlinear Lp - Lfir relation with a slope < 1 . For interact- 
ing g alaxies dSolomon & Sagd fl 9881) and LIGs dSanders et alj 
1991) the slope is approaching unity. In this case a linear rela- 
tion between CO and FIR emission can be assumed. For ULIGs 
ISolomon et alj d 19971) also find a linear relation between FIR and 
CO emission, with a higher Lfir/Lco ratio than in normal galax- 
ies and LIGs. The incr ease of the Lfir/Lco r atio towards ULIGs 
is well known (e.g.. [Sanders et al] Il986t ISanders & Mirabel 
1996: iGao & S olomon 2004). I t exte n ds towards the even more 
luminous SMGs dGreveetd] 120051) . IGao & Solomon] d2004l) 
point out, that a linear relation between CO and FIR luminos- 
ity remains valid only over two orders of magnitude in the lower 
L F ir range (while a corresponding HCN - FIR relation extends 
over three orders of magnitude and includes the ULIG domain). 
Any fit to th e CO luminosity over a wider range naturally yields 
a slope > l. lGao & Solomonl st ate a slope of 1.25 ±0.08 for L F1R 
between 10 95 L o and 10 12 5 L n . lYao et all d2003l) even obtain a 
slope of 1 .7 for Lfir between 10 9 L G and 1O 12 L O .0 

2 Their small projected beam size, however, may not cover the full 
CO content especially in the nearby, less infrared luminous galaxies. 
Thus, the slope may appear steeper than it actually is. 



In interferome tric CO data on nearby ULIGs 

iDownes & Solomonl d!998l) show that most of the emis- 
sion originates from rotating circumnuclear disks or rings. 
The disk gas forms a continuous medium rather than discrete 
virialized clouds. The linearity between FIR and CO emission 
in ULIGs and the dense environment resulting from the massive 
concentration of molecular gas suggest thermal excitation 
and optically thick emissio n. This is proposed by the same 
authors in an earlier paper ([Solomonet al.l 1 1 997b but partially 
revised in IDownes & Solomon (1998J): According to their 
model fits the CO line emission is only moderately opaque 
and the CO(2-l )/CO( 1-0) ratio indicates sub-thermal excitation. 

If a connection between ULIGs and QSOs exists, as sug- 
gested in the evolutionary model, what does this imply for the 
molecular gas content in the evolved QSO phase of these ob- 
jects? The model assumes a dispersion of the majority of gas 
and interstellar dust to clear the view onto the formerly dust en- 
shrouded AGN. Indeed, low-z PG QSOs and nearby HE QSOs 
show CO and FIR luminosities and star formation efficiencies 
that are lower than it is the case for ULIGs in the same comov- 
ing volume. Nevertheless, the majority of nearby low-luminosity 
QSOs is associated with large reservoirs of molecular gas (Sect. 
15.11 ) and presumably also with ongoing star formation, as it is 
indicated by their L' CQ - Lfir distribution (Fig. [8j. This might 
imply that the density of the formerly AGN obscuring medium 
in the central region is significantly reduced. 

But how does the evolution affect the distribution of molec- 
ular gas? A confinement of the molecular gas to a compact re- 
gion i s still indicated also in high redshift QSOs dRiechers et al.l 
2006, in ground-state transition observations at z^4). Spatially 
resolved data on QSOs i n the local universe , howe ver, is hardly 
available. For Mrk 23 1 IDownes & Solomonl d!998l) report on an 
inner and an outer face-on disk within a radius of 1.15 kpc. 
Stag uhn et alj d2004l) present a ringlike structure with a radius 
of ~ 1.2 kpc in the circumnuclear molecular gas distribution of 
IZwl . For the low-luminosity QSO HE 1029-1 83 l lKrips et all 
(2007) estimates a size of 6±2 kpc for the CO source that is 
aligned with th e optical bar. The in terfe rometric maps of the l ow- 
z PG QSOs bv lEvans et all d200ll) and lScoville etail d2003l) do 
not permit to draw conclusions on the compactness of the emis- 
sion. 

The coincidence with the LIG population in the logL c „- 
log L F ir diagram with a power-law index 1 can be taken as indi- 
cator for a confinement of remaining H2 content to a small and 
dense region. Based on this assumption, we would like to revisit 
the idea of both bands being optically thick. If this is the case, 
the two different inclinations in the log L co -log Lfir diagram 
may have a simple geometrical explanation: For normal galax- 
ies, the CO distribution is less confined. The totality of individ- 
ual gas clouds with varying line-of-sight velocities contributes 
to the spectrum. In this case, it is appropriate to use the velocity- 
integrated CO line emission as tracer for the total molecular gas 
content of the galaxy. The FIR luminosity, on the other hand, 
is dominated by reradiated dust emission, which originates from 
the densest regions in the ISM, the seeds of star formation. If 
the column densities in these compact regions result in t > 1 at 
the FIR bands, only the radiating surface contributes to the de- 
tectable FIR flux density. Let us assume a coarse, linear relation 
between the molecular gas content and the number of star for- 
mation seeds in a 3-dimensional gas distribution. Any surface of 
optically thick emission then scales with power of | of the over- 
all volume. If L' co can be used as measure for the total molecular 
gas content, Lfir should scale with L^ 67 and produce the corre- 
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sponding inclination in the log-log plot. Once the density of the 
ISM reaches a threshold at which the CO emitting region (which 
may consist of a continuous medium) becomes opaque, both the 
CO and the FIR emitting surface scale linearly with the overall 
molecular gas and dust content. Under these circumstances the 
slope in the log-log plot is close to unity. An even higher density 
of star formation seeds in the ISM results in a larger surface fill- 
ing factor within the galaxy, allowing for a higher Lpm/L' co ratio 
for each given CO luminosity. This effect causes the dispersion 
perpendicular to the power-law fit. As shown in the Appendix 
the different power-law fits for the individual source classes can 
be parameterized via the fractional source filling factor a of the 
material that participates in intense star formation or the source 
type / (QSO hosts and LIGs/ULIGs or normal galaxies). 

The LIG/QSO scenario and the normal galaxy scenario de- 
scribe the two extreme cases with power-laws 1 and 0.67 as up- 
per and lower limit. These extreme cases seem to correspond 
well to the two populations described above. Furthermore, these 
scenarios do not exclude a transition region, in which the dense 
central region of a galaxy is optically thick for both bands while 
for the outer regions the full volume contributes to the detected 
CO emission. 

6. Summary 

This paper reports on 12 CO observations of 41 nearby low- 
luminosity QSOs. Thirty-nine of them form a subsample that 
is not biased towards FIR emission. The results of this study can 
be summarized as follows: 

1. 70% of the subsample have been detected in the 12 CO(1-0) 
transition. With a mean H2 mass of 2.8 • 1O 9 M (assuming 
a — 4 M Q (K km s _1 pc 2 )~') the majority of low-luminosity 
QSO hosts is rich in mo lecular gas. This co nfirms previous 
results on PG QSOs bv lScoville eTail d2003h . 

2. For the majority of detected objects, the redshift based on 
visible features is blueshifted with respect to the CO line 
centroid - probably as a result of asymmetric phenomena in 
the vicinity of the nucleus, the emission region of the visible 
features. 

3. The absolute Bj magnitude of the central seeing disk of the 
low-luminosity QSOs is used as an indicator for the total 
activity of AGN and circumnuclear starburst. Objects with 
M(B j) < -20mag not only show higher CO luminosities but 
also a wider spread in the distribution of linewidths than the 
fainter objects. 

4. The connection between CO and FIR properties is discussed 
in the context of normal galaxies, local Seyfert I objects, 
LIGs, ULIGs and PG QSOs. Two populations with different 
power-laws in the L' cq -Lpir plot are identified in the range 

between Lfir = 10 9 L o and Lfir = 5 ■ 1O U L . The first 
population shows a power-law index of ~ 0.7 and contains 
normal, non- or weakly interacting galaxies and the bulk of 
Seyfert I objects. The second population with a power-law 
index of ~ 1 contains LIGs and the PG QSOs. The low- 
luminosity QSO sample seems to consist of objects of both 
populations and objects from a transition region. 

5. For Objects of the second population, i.e. also for galaxies 
hosting brighter QSOs, the linear relation between FIR and 
CO luminosity can be taken as indication for a concentration 
of the molecular gas in compact regions, similar to the case 
of ULIGs. 

6. The idea of thermalized, optically thick FIR and CO emis- 
sion, previously suggested in the context of ULIGs, is revis- 



ited. The two different slopes in the L^ -Lfir plot are dis- 
cussed as a potential result of different optical properties of 
the ISM in the two populations. 

Acknowledgements. TB would like to thank the staff of the IRAM 30m tele- 
scope for their kind support. We thank the anonymous referee for helpful com- 
ments. Part of this work was supported by the Deutsche Forschungsgemeinschaft 
(DFG), project number: SFB 494. 



Appendix A: Descriptive star formation model 

Acknowledging the opaqueness of a dense, potentially continu- 
ous, medium and the importance of the contribution of individ- 
ual star formation regions to the overall emission leads to the 
consequence that the sources are characterized by the surface 
filling factor a of the star forming molecular gas in the host and 
a quantity B that describes the Lfir and L' co contributions of the 
emitting regions (see below). Here L' CQ alone is not necessar- 
ily an adequate measure of the total molecular gas content in 
dense environments (as outlined in Sect. 15.4b . Only if the stan- 
dard assumptions are valid (virialized, i.e. gravitationally bound 
giant molecular clouds with different line-of-sight velocities), 
Ico traces molecular gas in both the LIG/ULIG/PG QSO case 
and the normal galaxy case. 

The two extreme populations (LIG/ULIGs/PG QSOs and 
normal galaxies) in Fig.|8]can then be described as follows. For 
PG QSO hosts and LIGs CO and FIR luminosity are directly 
proportional: 

lo g( L co) a log(L FIR ). 

For normal galaxies, additional quiescent molecular gas reser- 
voirs contribute to the object's total CO luminosity but do not, 
or at least not as efficiently, produce dust emission in the FIR. 
The additional contribution grows non-linear with increasing 
total molecular mass. The portion of CO-luminosity per FIR- 
luminosity is higher by a factor (1 +B), with B > 0: 

log(L^ ) cx [1 + y 0]log(L FIR ) ~ 1.41og(L FIR ). 

As a first step we assume in the following a smooth transition 
between the host properties of these different classes, which can 
be parameterized by the value / with 0</<l. f-l represents 
the case of QSO hosts and LIGs/ULIGs, while /=0 represents 
the case of normal galaxies. Furthermore, a direct proportional- 
ity between the CO and FIR luminosity in the /=1 case and in 
the star forming portions of normal galaxies is assumed. We also 
assume that in the f- 1 case the molecular gas of a dominant part 
of the host participates in the star formation process. In normal 




a~1.0 1.0<a<0.1 a~0.1 



Fig. A.l. Examples of molecular gas distributions with different 
fractions of molecular gas taking part in active star formation 
(dark areas). We show an inclined host and the outer continuous 
line includes the majority of the molecular gas in the system. 



T. Bertram et al.: Molecular gas in nearby low-luminosity QSO host galaxies 



13 




"i r 

0.5 0.8 1.0 
source type / 



Fig. A.2. Logarithm of the surface filling factor a of star-forming 
gas as a function of source type / as described in the text. The 
black filled circles mark the fits to the LIG/ULIG and the normal 
galaxies, respectively. The white filled circle indicate the region 
with which the nearby low-luminosity QSO hosts (discussed in 
this paper) are located. 



galaxies only a fraction a - for which log(L^ ) cc log(L F i R ) is 
valid - participates in that process. The quantity a can be identi- 
fied as the surface filling factor of the star forming molecular gas 
in the host with 0<a<l (cf. Fig. IA. 11 1. The remaining part gives 
rise to j81og(L F i R ) as an additional portion of CO-luminosity. 
With a parameter Co we can then write: 

log(Lco) - /log(L F iR) + 

(1 - /)[log(aL FIR ) +/81og((l - a) Lhr)] + C 

= [l+(l-/) y 6]log(L FIR ) + 

(l-/)[log( fl ) +y 01og(l-fl)]+Co 



In Fig. IA.2l we show as a second step the interdependency of a 
and / as they result in suitable power-law fits to the extreme and 
intermediate source classes in Fig. [8] We find that both quantities 
are linked via 

log(a) = 3(f-l). 

Within this model the two extreme source populations are rep- 
resented by the following cases: 

Case I: for f=l, i.e. for QSO hosts and LIGs/ULIGs with a^l 
follows: 

lo g( L Co) - l°g(^FIR) + Co . 

With a value of Co — 2, as suggested by Fig. [8] this corresponds 
to the heuristically derived expressions for LIGs/ULIGs and PG 
QSOs dSanders et alJfT99ll : ISolomon et dill 9971: IScoville et ail 
l200HlEvans eTaIll2006h . 

Case II: For normal galaxies with /=0, i.e. 0<a<Kl, and /3~0.4 
the model yields: 

log(L')s 1.41og(L FIR )-5.0. 



This corresponds to the heuristically derived expressions for 
samples of distanc e limited or non-interacting, isol ated FIR- 
selected galaxies dSolomon & Sagdll988t ISagdl 1993L cf. Sect. 
El. 

We see that the formalism described above naturally results 
in the heuristic relations for the extreme source populations as 
presented in section 5.4 and shown in Fig. [8] It also allows us to 
describe the mean properties of nearby QSO hosts discussed in 
this paper as an intermediate population. The intermediate case 
of the nearby QSO hosts, for which 0.5^/^0.8 or log a — 1, can 
then be approximated by 

log(Z^ ) a (1.1 ± 0.05)log(L FIR ) - (2.5 ± 0.3) 



or 



L FIR /10 y = (50±20)-[L'/10 y ] 



,9-1 0.91 (±0.04) 



Overall the model corresponds to a picture in which all hosts 
are similar in the amount of molecular gas (i.e. about 10% of 
~ 10 12 M Q of the total host mass) but different in the fraction of 
strong star forming regions. It also shows that the surface filling 
factor a and the value /3 that is linked to the FIR/CO emission 
properties of the source components are suitable parameters to 
describe a broad range of host galaxies. 
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